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Abstract 

We construct globally regular gravitating Skyrmions, which possess only discrete 
symmetries. In particular, we present tetrahedral and cubic Skyrmions. The SU(2) 
Skyrme field is parametrized by an improved harmonic map ansatz. Consistency then 
requires also a restricted ansatz for the metric. The numerical solutions obtained 
within this approximation are compared to those obtained in dilaton gravity. 



1 Introduction 

Nonlinear field theories coupled to gravity lead to globally regular gravitating configura- 
tions [1]. Moreover, also black hole solutions with nonlinear hair arise [1]. These black 
hole solutions are asymptotically flat and possess a regular event horizon. Outside their 
horizon they retain the features of the corresponding gravitating solitons, and may thus 
be viewed as bound states of solitons and Schwarzschild black holes [2]. 

In the Einstein-Skyrme model the nonlinear chiral field theory describing baryons and 
nuclei in terms of solitons (so-called Skyrmions) is coupled to gravity. Static spherically 
symmetric SU(2) gravitating Skyrmions and black holes with Skyrmion hair [3]- [6] exhibit 
a characteristic dependence on the coupling parameter: two branches of solutions merge 
and end at a maximal value of the coupling parameter, and only the solutions on the lower 



* Email: 
''Email: 

* Email: 



ti3@auth.gr 

klcihaus@theorie.physik.uni-oldenburg.de 
kunz@theorie.physik.uni-oldenburg.de 



branch are classically stable. Spherically symmetric SU(N) gravitating Skyrmions and 
black holes have been analyzed similarly [7, 8]. 

Recently, static axially symmetric SU(2) gravitating Skyrmions and black holes with 
Skyrmion hair have been constructed numerically [9], while approximations to axially sym- 
metric gravitating SU(3) Skyrmions have been obtained in [10]. The hairy black holes rep- 
resent further examples demonstrating that Israel's theorem does not generalize to theories 
with non-Abelian fields [11]. 

On the other hand, in flat space also Skyrmions with no rotational but only discrete 
symmetries have been constructed [12]. Among them are solutions with the symmetries of 
the platonic solids, to which we refer as platonic Skyrmions. Besides the exact numerical 
solutions, also approximate solutions have been obtained. Such approximate Skyrmion 
solutions with baryon number B have been constructed using rational maps of degree B 
between Riemann spheres [13], as well as an improved harmonic map ansatz [14]. 

In this paper we consider gravitating Skyrmion configurations with only platonic sym- 
metries. In particular, we focus on configurations with tetrahedral and cubic symmetry, 
possessing baryon number B = 3 and B = 4, respectively. 

Recall that the SU (2) Einstein-Skyrme action reads 



where R is the curvature scalar, g denotes the determinant of the metric, the 577(2) Skyrme 
field U enters via = d^UU^ 1 , G represents Newton's constant, and k and e are the 
Skyrme model coupling constants. 

In order for finite-energy configurations to exist the Skyrme field must tend to a constant 
matrix at spatial infinity: U — > 1 as — > oo. This effectively compactifies the three- 
dimensional Euclidean space into S 3 and implies that the Skyrme fields can be considered 
as maps from S* 3 into SU(2). 

Variation of action (1) with respect to the metric g^ v leads to the Einstein equations 
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with stress-energy tensor 
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The Einstein-Skyrme system has a topological current 
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which corresponds to the baryon current, yielding the baryon number B 
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For comparison, we also consider Skyrme-dilaton theory with action 
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with 7 = ^/AnG, and indices are lowered and raised by the flat Minkowski metric. 

The idea of the harmonic map ansatz for SU(N) Skyrmions [15] (which is the gener- 
alisation of the rational map ansatz of Houghton et al. [13]) involves the separation of the 
radial and angular dependence of the fields as 
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Here h(r) is the corresponding profile function and P is a N x N hermitian projector, 
which depends only on the angular variables (£,£), where £ is the Riemann sphere variable 
given by £ = e tlp tan(0 / '2) in terms of the usual spherical coordinates r, 9, (p. Note that the 
matrix P can be thought of as a mapping from S 2 into CP 1 * -1 . Thus, P can be written as 



P(V) = 



7^ 
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where V is a iV-component complex vector (depending on £). 

The improved harmonic map ansatz is obtained by allowing the profile function h to 
depend on all spatial coordinates r, As has been shown in [14], this ansatz leads to a 
better approximation of the Skyrmion energy in flat space. 

For SU(2) Skyrmions the harmonic map ansatz can be related to rational maps 

by 

V = (i2(0,l)*, (10) 

where = £>(£)/?(£) > an< ^ p(0> ?(0 are polynomials of £. Defining the unit vector 
by 

(11) 



(R + R, -i(R-R) , 1 - \R\ 2 ) , 



i + \R. 

the Skyrme field C/ takes the simple form 

U = cos(h) 1 + i sm(h)nR ■ f , 
with the vector of Pauli matrices f = {t\ , r 2 , r 3 ). 



(12) 
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2 Metric Ansatz 



To obtain static solutions without rotational symmetries let us first consider the following 
ansatz for the metric 

ds 2 = -fdt 2 + j (m x dr 2 + m 2 r 2 d6 2 + lr 2 sin 2 9dip 2 ) , (13) 

where /, mi, m 2 and I are functions of r, 9 and ip. In the case of axial symmetry, m x = 
m2 = m [16]. 

Insertion of the Skyrme field (12) and the metric (13) in the action (1), and subsequent 
variation of the action with respect to the Skyrme profile function h yields a second order 
partial differential equation (PDE) for h. Similarly, PDEs for the metric functions are 
obtained, which are equivalent to those obtained from the general Einstein equations (2) 
after insertion of the Skyrme field and the metric. 

It now turns out, that this coupled system of PDEs does not possess a solution in 
general. While surprising at first, this fact has the following reason: The Einstein tensor 
G^y is defined so that its covariant divergence vanishes, 

V M G^ = . (14) 

Consequently, the Einstein equations require 

V M T^ = , (15) 

provided there exists a solution. Insertion of the Skyrme field (12) and the metric (13) 
yields, however, 

V M T^° = , V M T^ r = , 

V M T^ ^ , V M T W ^ . (16) 

Consequently, not all Einstein equations can be satisfied. In the axially symmetric case, 
when mi = m 2 and the Skyrme and metric functions do not depend on tp, we can achieve 
V /j r w = 0, but still V^T^ 61 7^ 0. Only in the spherically symmetric case, when mi = m 2 = 
I and the Skyrme and metric functions depend only on r, all equations in (15) are satisfied. 
In this case, however, the ansatz for the Skyrmion field leads to an exact solution. 

This problem can be traced back to the fact that the harmonic map ansatz for the 
Skyrme field is too restrictive. Indeed, for a general ansatz for the Skyrme field U involving 
three functions hi(r, 9, ip) 

U = cos(/ii)l +ism(hi) [sin(/i 2 )(cos(/i3)Ti + sin(/i 3 )r 2 ) + cos(h 2 )r 3 ] (17) 

all equations in (15) are satisfied. 

While aiming at the numerical construction of exact platonic gravitating Skyrmions, in 
view of the complexity of the coupled Einstein-Skyrme equations we first want to obtain 
simpler approximate solutions, based on the improved harmonic map ansatz. 
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We therefore argue as follows: since we employ only an approximate ansatz for the 
Skyrme field, we should restrict also to an approximate ansatz for the metric, compatible 
with the ansatz for the Skyrme field. An appropriate ansatz for the metric is given by 

d g* = _ fdf 2 + \_ ^2 + r 2 M 2 + r 2 gin 2 ? (lg) 

where we allow only for two metric functions / and I (as in the spherically symmetric 
case), which, however, depend on all three coordinates, like the Skyrme function h of the 
improved harmonic map ansatz. 

We then derive the set of three coupled PDEs as variational equations from the Einstein- 
Skyrme action (1), after the second order derivatives of the metric functions have been elim- 
inated by intergration by parts. We refer to this approximation as "/ — /-approximation" . 

We note that a further restriction of the metric obtained by setting 1 = 1, leads to 
solutions of the Skyrme-dilaton model with action (7). In this case the dilaton can be 
expressed by the metric function /, 

= --Llog(/), 7 = v / 4^G (19) 

We refer to this approximation as "dilaton-approximation" . 

To obtain an estimate for the quality of the approximation, we can substitute the 
approximate solutions in the full set of Einstein and matter equations, to see how strongly 
these are violated. 



3 Numerical Solutions 

3.1 Parameters and Boundary Conditions 

Introducing the dimensionless radial coordinate x = Ker and coupling parameter a 
AttGk 2 , the action (1) becomes 



s = 5 / + i Tr {K » Kn + h Tr {[K ^ KA [K "> RU]) 



-gd 4 x, (20) 



and the Einstein equations read G^ v = 2aT fMU . Thus the solutions depend only on the 
coupling parameter a and the chosen rational map. Likewise, in Skyrme-dilaton theory, 
the solutions depend only on the dimensionless coupling parameter a = rfn 2 and the 
chosen rational map. 

For Skyrmions with axial symmetry rational maps of degree B are simply given by 
Rb = C, B , while for Skyrmions with platonic symmetries these maps are more complicated. 
In particular, for tetrahedral B = 3 Skyrmions and cubic B = 4 Skyrmions these maps are 
given by 

tftetra = ^ ( 21 ) 
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and 

e + gvgg + 1 
cubc " - 2v/3« 2 + 1 ' 1 j 

respectively. 

In order to map the infinite range of the radial variable x to the finite interval [0, 1] we 
introduce the compactified radial variable x — x/(l + x). 

For axially symmetric solutions the Skyrme and metric functions depend only on the 
coordinates x and 9. Due to the reflection symmetry, z <-> —z, it is sufficient to construct 
solutions for < x < 1, < 9 < n/2. The boundary conditions at the origin are 

h(0) = n, d x f\ = 0, d x l\ o = 0. (23) 

Asymptotically the Skyrme field tends to the unit matrix and the metric approaches the 
Minkowski metric, i.e. 

h(oo) -> , /(oo) -> I , Z(oo) -> 1 . (24) 

On the z-axis (9 = 0) and in the xy-plane (0 = 7r/2) the boundary conditions follow from 
regularity and reflection symmetry, respectively, 

d e h = 0, d e f = 0, d e l = 0. (25) 

For solutions with discrete symmetries the Skyrme and metric functions depend on 
all three coordinates x, 9, tp. The tetrahedral symmetry of the B = 3 solution allows to 
restrict to the domain of integration to < X < 1, < < 7r, 0<y?< tt/2. Similarly, the 
cubic symmetry of the B = 4 solution allows to restrict to the domain of integration to 
< x < 1, O<0<7r/2, < (p < ir/2. The boundary conditions at the origin, at infinity, 
on the z-axis and in the xy-plane are the same as for the axially symmetric solutions. The 
remaining boundary conditions at <p — and ip = n/2 follow from the platonic symmetries, 
i.e. 

d^h = , dj = 0, <y = . (26) 



3.2 Numerical Results 

Solutions are constructed with help of the software package FIDISOL [17] based on the 
Newton-Raphson algorithm. Typical grids contain 70 x 30 points for the axially symmetric 
solutions and 70 x 25 x 25 points for the platonic solutions. The estimated relative errors 
are approximately 0.1%, except close to a max , where they become as large as 1%. 

We have constructed gravitating Skyrmions with baryon number B = 2,3,4 in the 
"/ — /-approximation" and "dilaton-approximation" , and studied their dependence on the 
coupling parameter a. The a-dependence of the axially symmetric (B = 2,3,4) and pla- 
tonic (B = 3, 4) Skyrmions is completely analogous to the a-dependence of the spherically 
symmetric B = 1 Skyrmions [5]- [7]. 
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^max 


rat. map 


"f-1" 


"dilaton" 




0.0318 


0.0294 




0.0266 


0.0246 


Rtetra 


0.0267 


0.0248 


i?4 


0.0231 


0.0214 


Rcube 


0.0234 


0.0218 



Table 1: The maximal value of the coupling parameter, a max , for axial (i?2, -R3, -R4) and 
platonic (Rtetra, Rcube) Skyrmions. 



Gravitating Skyrmions exist only up to a maximal value of the coupling parameter, 
« max , which depends on the specific rational map and on the approximation (see Table 
1). When a is increased from zero a branch of gravitating Skyrmions emerges from the 
corresponding flat space Skyrmion solution. This first (lower) branch extends up to the 
maximal value a max , where it merges with a second (upper) branch of solutions. The upper 
branch then extends back to a = 0. 

In Fig. 1 we present the mass per baryon number as a function of a (left), for axially 
symmetric B = 2 and platonic B = 3 and B = 4 Skyrmions. On both branches the mass 
decreases with increasing a. But whereas the mass remains finite in the limit a — ■> on 
the lower branch, it diverges in this limit on the upper branch. Thus on the upper branch 
the limit a — > does not correspond to a flat space limit, where gravity decouples. 

To better understand the limit a — > on the upper branch, we note that the coupling 
parameter a = AtiGk 2 vanishes, either when G vanishes while k remains constant, or when 
k vanishes while G remains constant. The first case corresponds to the flat space limit of 
the lower branch, while the second case corresponds to the limit a — > of the upper branch. 
Introducing the rescaled radial coordinate x = x/y/a and the rescaled mass M = My/a, 
one observes, that the rescaled mass remains finite in the limit a — > on the upper branch, 
as illustrated in Fig. 1 (right). 

We note that the a-dependence of the mass of the gravitating Skyrmions solutions is 
almost the same in the "/ — /-approximation" as in the "dilaton-approximation" . In the 
"/ — /-approximation" the mass is slightly higher, in particular, along the upper branch. 
Also, a max is slightly larger in the "/—/-approximation" . The mass of the axially symmetric 
B = 3 and B = 4 Skyrmions is always larger than the mass of the corresponding platonic 
Skyrmions. 

In Fig. 2 we exhibit the value of the metric functions / and / at the origin for these 
axial (B = 2) and platonic (B = 3, 4) Skyrmions. We observe that /(0) and 1(0) take finite 
values in the limit a — > on the upper branch. 

Focussing now on the platonic Skyrmions, we exhibit in Figs. 3 surfaces of constant 
baryon density for tetrahedral B = 3 and cubic B = 4 Skyrmion solutions on the lower 
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0.04 



Figure 1: The dimensionless mass per baryon number M/B (left) and the scaled mass per 
baryon number (M/ B)*Ja (right) are shown as functions of the coupling parameter a for 
axial (B = 2) and platonic (B = 3, 4) Skyrmions in the "/ — /-approximation" and the 
"dilaton-approximation" . 

branch (left) and upper branch (right). For a given rational map and coupling parameter 
a the Skyrmion on the upper branch is confined in a smaller volume than the Skyrmion 
on the lower branch, while the shape of the baryon density is primarily determined by the 
rational map [12], analogous to the shape of the energy density [18]. 

In Figs. 4 we demonstrate that the metric functions / and / of the platonic B = 3 and 
B = 4 Skyrmions possess the same symmetry as the baryon density. In fact, when gravity 
is weak the function 1 — / is proportional to the Newtonian gravitational potential. 




0.04 



Figure 2: Same as 1 for the value of the functions / and / at the origin. Note, that I = 1 
in the "dilaton-approximation" . 
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4 Conclusions 



While aiming at the numerical construction of exact platonic gravitating Skyrmions, we 
have here obtained simpler approximate solutions, based on the improved harmonic map 
ansatz for the Skyrme field, thus avoiding the full complexity of the coupled Einstein- 
Skyrme equations. This ansatz for the Skyrme field involves a single function instead 
of three functions. Consequently, an appropriate restriction of the ansatz for the metric 
is required, involving either two functions ("/ — /-approximation") or a single function 
( "dilaton-approximation" ) . 

We have focussed on platonic gravitating Skyrmions with tetrahedral (B = 3) and 
cubic (B = 4) symmetry. For comparison, we have also constructed gravitating axially 
symmetric Skyrmions (B = 2,3,4). The dependence on the coupling parameter a of the 
axial and platonic Skyrmions is completely analogous to the a-dependence of spherical 
B = 1 Skyrmions. When a is increased from zero a lower branch of gravitating Skyrmions 
emerges from the corresponding flat space Skyrmion solution. This branch extends up to a 
maximal value a max , where it merges with an upper branch of solutions, which extends back 
to a = 0. Thus gravitating Skyrmions exist only up to a maximal value of the coupling 
parameter, a max , which depends on the specific rational map and on the approximation 
used. 

The shape of the baryon density of platonic Skyrmions is primarily determined by the 
rational map and analogous to the shape of the energy density. For a given rational map 
and coupling parameter a the Skyrmion on the upper branch is confined in a smaller volume 
than the Skyrmion on the lower branch. The metric functions of the platonic Skyrmions 
possess the same symmetry as the baryon density. 

Comparing the approximations, applied in the construction of the gravitating tetra- 
hedral and cubic Skyrmions, we note that their mass is slightly higher in the "/ — l- 
approximation" than in the "dilaton-approximation" . On the other hand, in flat space the 
exact Skyrmion [12] has a slightly lower mass than the approximate Skyrmion [14]. There- 
fore, it is an open question whether the mass of the exact gravitating Skyrmion solutions 
remains lower than the mass of the approximate solutions for all values of the coupling 
parameter. Construction of the exact platonic gravitating Skyrmion solutions, however, 
remains currently still a numerical challenge. 
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Figure 3: Isosurfaceplot of the baryon density B°(x, y, z) = 0.25 x max(_B°) for the B = 3 
Skyrmion (upper row) and B°(x,y,z) = 0.5 x max(i? ) for the B = 4 Skyrmion (lower 
row) in the "/ — /-approximation" for a = 0.02 on the lower branch (left column) and the 
upper branch (right column). 
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Figure 4: Isosurfaceplots of the functions 1 — / (left column) and 1 — / (right column) 
for the B = 3 Skyrmion (upper row) and the B = 4 Skyrmion (lower row) in the "/ — l- 
approximation" for a = 0.02 on the lower branch. 
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